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The Sun is a typical main sequence star that generates its energy by the pp chain 
and CNO nuclear cycle , which also produce neutrinos. These neutrinos have been 
observed by four pioneering solar neutrino (uq) experiments, the Chlorine exper- 
iment at Homestake, the Cerenkov light water experiment at Kamiokande, the 
Soviet-American Gallium experiment (SAGE) at Baksan and the European Gal- 
lium experiment (GALLEX) at Grand Sasso, but the measured uq fluxes are signifi- 

cantly below those predicted by Standard Solar Models (SSM): The CI experiment 
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observed over 20 years an average production rate of 2.28 ± 0.23 SNU {SNU = 
10^"^^ captures/atom s) of "^"^Ar by solar u^s with energies above 0.81 MeV while 
the SSM predicts 7.4 ± 2.8 SNU (unless otherwise stated, we will use the predic- 
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tions of Bahcall and Pinsennault (BP) as a standard reference). Kamiokande has 
observed since 1987 electron recoils with energy above 7.5 MeV from scattering 
of z/q's during 1040 (Kamiokande II) and 395 (Kamiokande III) days of detector 
live time with event rates that are 51% ± 6% and 56% ± 6%, respectively, of that 
predicted by the SSM. GALLEX has measured a capture rate of solar z/g's by ^^Ga 
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of 87± 16 SNU in 21 runs between May 1991 and February 1993, compared with 

128+1^ SNU, predicted by the SSM. SAGE began its uq observations in 1990. Its 
initial results indicated a large discrepancy with the SSM predictions . However, 
the capture rate observed by SAGE in 1991 , 85^32[stat] ± 20[si/s] SNU, agrees 
with that observed by GALLEX. In spite of impressive theoretical and experi- 
mental efforts the origin of these discrepancies, which became known as the uq 
problem, is still unknown. 

o 

Bahcall and Bethe have argued that the solution of this uq problem requires new 
physics beyond the Standard Electroweak Model because the CI detector with an 
energy threshold of 0.81 MeV observes a smaller solar z/g flux than that observed 



by Kamiokande whose energy threshold is 7.5 MeV. However, after the instalment 
of new pumps in 1986, the average capture rate measured by the CI detector is not 
in a serious conflict with the uq flux measured by Kamiokande during the same 
period (see Fig.l), especially in view of the absence of direct calibration of the CI 
experiment. Moreover, the uq spectrum measured by Kamiokande II and III is 
consistent with the z/g spectrum of ^B, as predicted by the SSM. Furtheremore, the 
observed capture rates in GALLEX and SAGE are those predicted directly from the 

solar luminosity and the observed fluxes of the energetic z/q's in the Homestake 
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and Kamiokande experiments, using essentially only general conservation laws : 

Since the net outcome in the pp and CNO cycles is the conversion of protons into 

Helium nuclei, conservation of baryon number, charge lepton flavour and energy 

requires that 4p + 2e^ -^ '^He + 2z/e + Q , where Q ~ 26.73 MeV, i.e., two z/g's 

are produced in the Sun per 26.73 MeV release of nuclear energy. Thus, if the Sun 

is approximately in a steady state where its nuclear energy production rate equals 

its luminosity then the uq flux at Earth is given by 

where Lq = 3.826(8) x 10"^"^ erg s^^ is the Sun's luminosity, D ^ 1.496 x 10^^ cm 
is its distance from Earth and Ei, ^ E,y{pp) ^ 0.265 MeV is the average z/q 
energy (Kamiokande and the CI experiments observe much smaller fluxes of the 
more energetic solar neutrinos). If the suppression of the flux of z/q's other than 
the pp neutrinos is energy independent and equals that observed at Homestake, 
S = (t)i^Q/<PsSM ~ 0.31, then the predicted capture rate in Galhum is < (pa >Ga~ 
77-65+575 SNU ^ 93 SNU, where 77-65 SNU is the capture rate of the pp z/q's 
and 575 SNU is the capture rate of all other z/q's. If the suppression is 5 ~ 0.51, 



as observed in Kamiokande, then < (pa >Ga^ 103 SNU. Both predictions are in 
good agreement with the observations of GALLEX and SAGE. All these suggest 
that perhaps the difference between the results of the CI experiment prior to 1986 
and those of Kamiokande after 1986 are due to statistical fluctuations or systematic 
errors, allowing a standard physics solution to the z/q problem. 

Indeed, in this letter we show that the SSM predictions are consistent with the z/q 
observations of Kamiokande, SAGE and GALLEX. This is achieved mainly by: (1) 
using an improved numerical stellar evolution code to calculate the evolution of 
the Sun from its pre-main sequence phase to the present stage, (2) using an updated 
solar luminosity Lq = 3.826(8) x 10^^ erg s~^ instead of 3.86 x 10"^"^ erg s^^ 
used before, which reduces the ^Be and ®B z/q fluxes by 3% and 6%, respectively, 
(3) using the updated values 6*34(0) = 0.45 keV ■ b and 5*17(0) = 17 eV ■ b for the 
low energy nuclear reactions ^He(a,7)'''Be and '''Be(p,7)^B, respectively, instead of 
'S'34(0) = 0.56 keV ■ b and 5*17(0) = 22.4 eV ■ b, (4) suppressing the overestimated 
screening enhancement of bare nuclear cross sections in the solar plasma (induces 
only a minor change after the above changes). 

The rate of the neutrino producing reactions in the Sun depends on the present 
conditions inside the Sun. These are obtained from calculating the evolution of the 
Sun from its pre-main sequence fully convective stage to its present (tQ ^ 4.55 GY) 
stage. Its initial composition, which is not known accurately enough, is adjusted 
by requiring that the calculations yield the presently observed properties of Sun, 

i.e., its radius, luminosity, observed surface elemental abundances and an internal 
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structure consistent with helioseismology data . The above scheme with the best 

available input physics (equations of state, opacities, nuclear cross sections) has 

been improved and updated continuously by Bahcall and his collaborators and 



more recently also by other groups . Our predictions below were obtained with an 
improved stellar evolution code which contains diffusion and partial ionization of 
all elements, follows the full evolution of the Sun from its pre-main sequence stage 
until present and does not impose nuclear equilibrium at any stage. The code uses 
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the thermonuclear reaction rates compiled by Caughlan and Fowler except for 

the changes explained below. 

Low energy nuclear cross sections in the SSM were extrapolated from mea- 
surements at much higher energies, using the parametrization, 

a{E) = ^e-(^-/^^''' (2) 

E 

where Eq = 2(naZiZ2)'^^c^, Z\ and Z^ being the charge numbers of the colliding 
nuclei, /x their reduced mass and E their center-of-mass energy. SiE) was assumed 
to vary slowly with energy . However, the above parametrization ignores the 
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non negligible contribution from partial waves other than s wave at the high 
energies where the cross sections have been measured and the effect of a finite 
nuclear radius R on barrier penetration, which induce an energy dependence in 
SiE). Indeed, SiE) was found to be energy dependent for some key reactions and 
its extrapolation to low energies became model dependent. However, the use of the 
correct Coulomb barrier penetration factor for finite size nuclei and the inclusion of 
higher partial waves remove most of the energy dependence of SiE). In particular, 
for effective potentials V = —Vq < for r < i? and V = ZiZ2e /r for r > R, the 
induced energy dependence of S{E) well below the Coulomb barrier and far from 
resonance energies is given by 



SiE) ~ S" " V ^ y ['2Vx-arcsin^~y/x{l~x)] 



where S = S{0) + S'{0)E, x = E/Ec and Ec = ZiZ2e^ / R . That is demonstrated in 

Fig. 2 for the reaction He(Q;, 7) Be, where we plotted S^4^{E) as function of E^ using 

R = i?(^He)+i?(a)=3.82 fm and the experimental data of Parker and Kavanagh 
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and of Osborne et al. Indeed, 6*34 shows very small energy dependence at low 
enough energies where the contribution from higher partial waves is negligible. 
With this form for the barrier penetration factor, the world experimental data on 
the reaction ^He(«, 7)^Be yield 6*34(0) = 0.45 ± 0.03 keV ■ b, which is significantly 
smaller than the value 6*34(0) = 0.56 keV ■ b used by Bahcall and Pinsennault . 
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Using the experimental data of Filippone et al. on the reaction ^Be(p,7)®B at low 
energies, where the contribution from higher partial waves is small, and the recent 
measurements by Motobayashi et al. of Si7{E) at higher energies from Coulomb 
dissociation of ^B, which are free of the resonating p-wave and higher partial waves 
that contribute to the law energy cross section of ^Be(p,7)®B, we find that for 
2.4 fm < R{p) + -R(^Be)< 3.2 fm, Sn is practically energy independent and 

5*17(0) = 17±2 eV -b (see Fig. 3). This value which we adopt was advocated before 
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by Barker and Spear and it is significantly smaller than 6*17(0) = 22.4 ± 2 eV ■ b 

2 
used by Bahcall and Pinsennault . 

Screening of the Coulomb potential of target nuclei by their electrons is known to 

enhance significantly laboratory nuclear cross sections at very low energies , al- 
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though a complete theoretical understanding of the effect is still lacking . Screen- 
ing corrections to the nuclear reaction rates in the SSM are usually represented by 
an enhancement factor 

p _ ^AU/kT _ ^Z,Z,e^/RnkT ^ (^^-^ 

where AU is the gain in electronic potential energy when an incident ion of charge 



Zje penetrates the electronic cloud of an ion of charge ZjC, T is the plasma tem- 
perature and Ry) is the Debye length, 



kT \ 



The change in the Coulomb barrier due to screening is estimated from the Debye- 
Huckel approximation to the screened potential around a static ion in an electrically 
neutral plasma, {TjZfiz = 0, nz being the number density of particles of charge Ze 
with Z=-l for electrons): 



$^ = £^e-^/^- ~ £^ _ ^ for r « Rd- (6) 

Rd r Rd 



This potential is an approximate solution to Poisson's equation near an ion of 
charge Zi, 

which is valid, however, only far from the nucleus where eZ^i <ti kT . In the 
core of the Sun where kT ~ 1 keV is much smaller than the Coulomb barriers 
(~ 1 MeV) between the reacting nuclei, most of the contribution to the nuclear 
reaction rates, comes from collisions with cm. energies E ^ kT . At the classical 
turning point eZj^i = E . Consequently, inside the barrier, eZj^i ^ kT, and 
the use of the Debye-Huckel approximate potential at distances shorter than the 
classical turning point for calculating the barrier penetration factor is unjustified. 
The Debye-Huckel solution also requires that the inter-ion spacing is much shorter 
than the Debye length, and that there are many electrons within a Debye sphere. 
Both conditions are not satisfied in the core of the Sun. It was found that when the 



conditions for the validity of the Debye-Huckel approximation are badly violated 
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in laboratory plasmas, the Debye-Huckel screened potential fails dramatically in 
reproducing various plasma properties. 

Moreover, an ion that approaches from infinity a nucleus screened by an electronic 
cloud does gain electronic potential energy At/ ^ e^ZiZj/ Ro when it penetrates 
this cloud. But, there is no gain in potential energy if either its initial position is 
already inside the electronic cloud (because the potential ~ eZi/ Rjj , is constant 
there), or the ion leaves and enters similar potential wells. Near the center of the 
Sun, where pc ^ 156 g cm-^, T^ ^ 1.57 xWK, Xc{H) ^ 34%, and X^iHe) ^ 64%, 
the average interion spacing is n^ ~ 2.8 x 10^^ cm, similar to the Debye length, 
Rd ~ 2.3 X 10^^ cm. Furthermore, the effective energies of the reacting ions 
are much higher than kT and their velocities are similar to the thermal electron 
velocities, not allowing sufficient time for the plasma to readjust itself and screen 
effectively these ions. Consequently, both effects reduce substantially the energy 
gain due to screening and the resulting enhancement of the fusion reaction rates 
near the center of the Sun. The overestimated enhancement of all nuclear cross 
sections by screening in the SSM has a small net effect on the Sun but may be 
important for other stars and will be further investigated in numerical simulations 
of stellar like plasmas . 

Our SSM predictions for the solar neutrino fluxes, obtained from the stellar evo- 
lution code of Kovetz and Shaviv with the above modifications are summarized 
in Table I. They are consistent, within experimental uncertainties, with the ex- 
perimental results of Kamiokande, GALLEX and SAGE, and of Homestake after 
1986. In particular, we predict 0,y(^B)=2. 77x10^ cm.^'^s^^ in excellent agreement 
with the fiux measured by Kamiokande II and III, but we do not explain why the 



uq flux measured by the CI experiment prior to 1986 is smaller than the uq flux 
measured later by Kamiokande. A full account of our calculations will be published 
elsewhere. 
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Table I: Comparison between experimental results and the predictions of the 
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SSM of Bahcall and Pinsennault (BP), Turck-Chieze and Lopes (LP), Kovetz 

and Shaviv (KS) and our predictions (DS) as described in the text, for the solar 

neutrino fluxes (cm^'^s^^) at Earth and their capture rates (SNU) in "^"^Cl and 

^^Ga. 

z/0 Flux BP TL KS DS EXP 

pp (ElO) 5.95 6.02 5.99 6.04 

pep (E8) 1.41 1.39 1.38 1.40 

^Be (E7) 4.75 4.34 4.91 4.30 

^B (E6) 5.24 4.63 5.83 2.77 2.87 ±0.17 

i^N (E7) 43.9 38.3 6.11 7.47 

1^0 (E7) 37.5 31.8 1.78 2.17 

i^F (E6) 4.72 4.85 5.21 

< 00- >ci 7.4 6.4 7.6 4.2 2.28 ± 0.23 

<(l)a>Ga 128 123 123 109 87 ±16 
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FIGURE CAPTIONS 

Fig.l: Comparison between the observed rates induced by solar neutrinos in the 
CI (open circles) and the Kamiokande II (square boxes) experiments divided by 
the rates predicted by the SSM of Bahcall and Ulrich (BU) as function of time 
during the period January 1987 through April 1990. The dotted line is the average 
of the CI data points, the full line is the average of the Kamiokande II data points. 

Fig. 2: 5*34 as function of Ecm as extracted from the cross section for ^He(Q;, 7)^Be 
measured by Parker and Kavanagh (triangles) and by Osborne et al. (full cir- 
cles). The straight line represents the average 6*34 obtained from the world data 
(not shown here). 

Fig. 3: 5^17 as function of Ecm as extracted from the low energy cross section for 
'''Be(p,7)®B measured by Filippone et al. (triangles) and from Coulomb dissocia- 
tion of ®B measured recently by Motobayashi et al. (full circles). 
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